Other incarnations of the Gross-Pitaevskii dark soliton 
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We show that the dark soliton of the Gross-Pitaevskii equation (GPE) that describes the Bose-Einstein 
condensate (BEC) density of a system of weakly repulsive bosons, also describes that of a system of strongly 
repulsive hard core bosons at half filling. As a consequence of this, the GPE soliton gets related to the magnetic 
soliton in an easy-plane ferromagnet, where it describes the square of the in-plane magnetization of the system. 
These relationships are shown to be useful in understanding various characteristics of solitons in these distinct 
many-body systems. 
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1. Introduction 

The study of solitary waves and solitons has been one of the most active areas of research in nonlinear systems for the past 
three decades (2. As is well known, a solitary wave is a spatially localized traveling wave which maintains its shape, size and 
speed as it travels. If a solitary wave retains all its features even after a collision with another solitary wave, it is called a soliton. 
Some of the iconic examples of completely integrable nonlinear partial differential equations arising in applied mathematics 
and exhibiting such strict soliton solutions O include the Korteweg-deVries equation, the sine Gordon equation, the nonlinear 
Schrodinger equation, etc. However, in physics literature a solitary wave itself is often referred to as a soliton due to its particle- 
like properties and we will do so in this paper as well. Solitons appear in diverse fields [1J such as hydrodynamics, fiber optics, 
biophysics, high energy physics, magnetism and in various condensed matter systems including a Bose-Einstein condensate 
(BEC). In addition to theoretical investigations, they have also been observed experimentally in various nonlinear systems. 

A BEC is an intrinsically nonlinear system. The Gross-Pitaevskii equation (GPE) (3) has been central to explaining various 
fundamental aspects of the BEC of a system of weakly repulsive bosons. One of the hallmarks of this system is the dark soliton 
that has been theoretically predicted as a unidirectional traveling wave solution for the condensate density. This soliton, 
which dies out when its speed approaches the speed of sound, has also been observed experimentally [5]. The study of solitary 
wave propagation in BEC remains an active frontier with emphasis on non-GPE dynamics [6| and many body effects such as 
quantum fluctuations and depletion that cannot be described by GPE (7J|8). 

In view of the relatively few examples of systems exhibiting solitons, any relationship that is found between different physical 
systems is fascinating, since it unveils the universal mathematical aspects underlying these nonlinear systems that describe 
completely different physical phenomena. In this paper, we show that the dark soliton of the GPE that describes the BEC density 
of a system of weakly repulsive bosons, also describes that of a system of strongly repulsive hard core bosons at half filling. This 
in turn leads to a connection between the GPE soliton and a magnetic soliton in an easy-plane ferromagnet. These relationships 
are shown to be useful in understanding various characteristics of solitons in these distinct many-body systems. 

The complex order parameter (i.e., the condensate wave function) describing a BEC that satisfies the GPE can be identified 
with the bosonic coherent state average [9] of the boson annihilation operator in the continuum version of the Bose-Hubbard 
model. Analogously, in a strongly interacting limit described by a system of hard-core bosons (HCB) the appropriate order 
parameter for this system would be the spin coherent state average [10] as the system maps to a pseudospin model [ 1 1 1. In the 
continuum description of the dynamics of the extended lattice Bose-Hubbard model for hard-core bosons with nearest neighbor 
interactions, this order parameter satisfies an evolution equation that is different from the GPE 1 1 2 1 . We will refer to this equation 
as the HGPE, where the prefix "H" has been used to denote its connection to HCB. 

Unlike the GPE which describes the dynamics of a BEC with no depletion, the HGPE evolution encodes both the normal and 
condensate fractions in the bose system. In our recent study, we showed that the HGPE supports 1 12 1 both a dark soliton and 
an antidark| 13 1 soliton (i.e., a bright soliton on a pedestal) for the bosonic density. In other words, the dark soliton of the BEC 
of weakly repulsive bosons acquires a partner that is bright, when the repulsion becomes extremely strong. Apart from being 
bright, the anti-dark soliton is found to be of quite a distinct variety compared to the dark soliton that resembles the GPE soliton. 
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Away from half-filling, i.e., when there is a nonvanishing particle-hole imbalance in the background, the dark soliton dies out 
when its propagation speed approaches the speed of sound (like the GPE soliton), while its brighter partner is found to persist 
all the way up to the sound velocity lfT2ll . 

In contrast, at half-filling, when the particle-hole imbalance is zero, the dark and the antidark solitons in the particle density 
profile of the HGPE become mirror images of each other, and they both die out at the sound velocity. Interestingly, in the 
corresponding condensate density, the two solitons become degenerate and cease to exist as separate entities. In this special 
case, the condensate density profile is a dark soliton that dies out as its speed approaches the speed of sound [ 12 1. 

In this paper, we show that for the half-filled case, the localized functional form of the solitary wave solution for the density 
of the hard-core bosonic condensate described by the HGPE agrees with that of the dark soliton of the GPE, to a very good 
approximation. In other words, if one focuses on the condensate fraction of the hard-core bosonic system, this strongly repulsive 
bose gas supports a solitary wave that has almost the same profile as the GPE soliton that describes the BEC of a weakly repulsive 
bose gas. 

Furthermore, since the HCB gas in the spin coherent state formulation of the extended Bose-Hubbard model with nearest 
neighbor interactions mimics an anisotropic, easy-plane ferromagnetic chain, the well known dark soliton of the GPE also relates 
to the magnetic soliton in the easy-plane ferromagnet, where it describes the profile of the square of the in-plane magnetization. 

In section 2, we begin with the many-body Bose-Hubbard Hamiltonian that describes the low-energy behavior of bosonic 
atoms in an optical lattice, and briefly outline some steps that lead to the order parameter equations in the weakly and strongly 
interacting limits, respectively. In sections 3 and 4, we first derive the evolution equation for the traveling waves in the 
condensate fraction of the atomic cloud obtained from the HGPE, and show that to a very good approximation, this equation 
agrees with the corresponding equation derived from the GPE. In section 5, we discuss the relationship between the density 
solitons of the BEC and the magnetic solitons of an easy-plane ferromagnet. 

2. The Bose-Hubbard model, the GPE, and the HGPE 

We begin with the extended lattice Bose-Hubbard model in d dimensions, whose Hamiltonian is given by 

H = -Y^[tb]b j+a + Vn,jn j+a ] + ^ U nj {n - 1) - (p - 2t)n h (1) 

3,0. j 

where and bj are the creation and annihilation operators for a boson at the lattice site j, nj is the number operator, a labels 
nearest-neighbor (nn) separation , t is the nn hopping parameter, U is the on-site repulsion strength, and [i is the chemical 
potential. To soften the effect of strong onsite repulsion, we add an attractive nn interaction {V > 0). Such a term may 
mimic certain characteristics of the long range dipole-dipole interaction, that has been considered in several recent studies [ 14] . 
The term 2trij is added so that the terms involving t reduce to the kinetic energy expression in the continuum version of the 
many-body bosonic Hamiltonian. 

Conventionally, one defines the order parameter for a many boson system to be the thermodynamic expectation value of the 
boson field operator. Invoking the concept of a broken gauge symmetry allows this expectation value to be nonzero below 
the BEC transition temperature. It has been argued|9| that this order parameter may be chosen to be the expectation value of 
the boson annihilation operator in the bosonic coherent state (also known as the harmonic oscillator coherent state or Glauber 
coherent state) representation of the pure quantum state. In this description, it can be shown that quantum fluctuations are absent. 

The Heisenberg equation of motion for the boson annihilation operator determined from ([lj, after a bosonic coherent state 
averaging yields the following GPE equation for the condensate order parameter ^(r, t) in the continuum description. 

- (fi 2 /2m)V 2 * 3 + U\V g \H g - = ihd t V g . (2) 

The GPE with U > provides a very useful characterization of the various properties of the BEC of weakly repulsive bosons 
101 . Here the condensate density p g = g \ 2 satisfies a continuity equation, and the system shows no depletion. 

The HCB limit (U — > oo) of the Bose-Hubbard Hamiltonian ([T]i has emerged as a useful model to describe various char- 
acteristics of the BEC of a strongly repulsive bose system. The constraint that two bosons cannot occupy the same site can 
be incorporated in the formulation by using field operators that anticommute at the same site but commute at different sites, 
thus satisfying the same algebra as that of a spin-^ system. By identifying bj with the spin flip operator Sj~, along with 
rij = Sj~Sj~ = | — Sj, the system can be mapped to the following quantum XXZ Hamiltonian in a magnetic field: 



H S = - J> Sj • S j+a - g S*S* +0 ] - - aO £/. 

3,0. j 



(3) 
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Here g = (t— V)d, where d is the spatial dimensionality. As we shall show, g = (t — V) > (see below Eq. (|8]l). Hence the 
HCB system maps to a quantum spin-1/2 system with a Heisenberg exchange interaction f, an easy-plane exchange anisotropy 
g and a transverse magnetic field h z = (g — //) along the z-direction. 

This mapping to spins suggests that a natural choice for the condensate order parameter ty s of the HCB system is the average 
of the spin flip operator in the spin coherent state representation lfT31l . i.e., 

* s =<S+>. (4) 

In this representation, it can be shown that lfl2l the condensate density, p s — | | 2 =< S~ >< S + > is related to total particle 
number density p =< S~ S + > as, 

Ps = P(l - P) (5) 

Setting ^> s = ^/pjexpi^ in the continuum description, the condensate order parameter ^> s and the particle density p satisfy 
the following equations l' T2l 

h 2 

iM s = - — {!- 2p)VH s - V e V s V 2 p + 2gp$ s - p$l s (6) 



p=r^V-[p(l-p)V$, (7) 
2m 

with the identification ta 2 = h 2 /m and Va 2 = V e . Equation (raj) (which can also be written as coupled equations for p and 0) 
will be referred to as HGPE. Using the asymptotic value p —> p in Eq. (|6J, the chemical potential p is found to be p — 2gp°. 

From Eq. Q, we note that it is important to distinguish between the total bosonic particle density p and the condensate density 
p s , where p = p s + pd- Here, pd describes the depletion, i.e., the normal component of p. 

Linearizion of HGPE by considering small amplitude fluctuations in the asymptotic value of ^ s leads to a Bogoliubov-like 
spectrum. For small momenta, the spectrum becomes linear, leading to the sound velocity c s for the hard-core system as 

c s = (2gp s /m)i, (8) 

where p° s is the asymptotic value of the condensate density. Since c s must be real, we obtain the condition 

g=(t-V)>0. 

In the next section, we will derive an equation of motion for the condensate density p s that arises from HGPE, for the 
half-filled case, when the asymptotic particle density is half, i.e., when the number of particles equals the number of holes in the 
background. 



3. HCB Condensate density evolution in the half-filled case of HGPE 



By seeking unidirectional traveling wave solutions (along the x-direction) of the total density p of the form, 

p(z)=p + f(z), (9) 

where z — (x — vt) j d and po denotes the constant asymptotic value of the background density, it has been shown iTPH that the 
HGPE (Eq. (|6|) leads to a nonlinear differential equation for j(z) which can be solved analytically to a very good approximation, 
to yield soliton solutions. 

Here we focus on the half-filled case with po = 1/2. For this case, the differential equation for / takes the form 1 12 1 



) 2 = 4/ 2 [ 7 2 -/ 2 ] 



(10) 



Here, 7 2 = 1 — v 2 with v — v/c s . Further, z — Qz and C, — A/Vl — 2A 2 , where the microscopic dimensionless parameter 
A = c s /cq is the speed of sound c s in the HCB condensate measured in units of the zero-point velocity cq — h/ma. 
As discussed in lfl2l . Eq. ( [T0| provides the following two soliton solutions 



f(z) = ±( 7 /2) sech(2 7 z) = ±( 7 /2) sech(z/r s ), 



(11) 



where the soliton width T s = (27C) 1 . Since 7 must be real in Eq. (Hi, the dimensionless parameter v must satisfy < v < 1. 
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The two solutions for / in ( [TT] ) lead to a doublet of localized solitons, namely, a dark soliton and an antidark soliton, for the 
density of the bosonic cloud p{z) = (1/2) + f(z). Using this in the condensate density expression given in Eq. |5]), we get 



Ps{z) 



(12) 



Writing p s (z) = (1/4) + f s {z), and comparing it with Eq. ( 12 1 shows that the variation f s around the asymptotic condensate 



density and the variation / around the total density p are related as follows: 



fs{z) = Ps{z) 



-f 2 (z) 



(13) 



Differentiating both sides of Eq. ( 13 i with respect to z and using Eq. ( |T0] > in the resulting equation, we obtain the following 
differential equations for the normalized condensate density p s = p s /(l/4) and the corresponding normalized variation f s = 
f s / (1/4), respectively. 



az 



(14) 



(^) 2 = i?(/ s +7 2 ) 



Unlike the doublet soliton solutions of / given in Eq. (Ill,/, has the unique solution 

f s = -4/ 2 = - 7 2 sech 2 2 7 z = - 7 2 sech 2 z/T s . 



(15) 



(16) 



Figure 1 shows the dark soliton profiles for the condensate density for various values of v. As the plot shows, the soliton 
which is dark when v = 0, dies out as v 1, i.e., as its propagation speed approaches the speed of sound. 

4. Condensate density evolution in the GPE 

In this section, we show that Eq. ( p"4| ), which describes the evolution of the condensate density of the HGPE for a half-filled 
system, has the same form as the condensate density evolution equation obtained from the GPE. 

Following the earlier studies of solitons in GPE|3 1, we write the condensate wave function in that system as 



' P°n (A + iv) 



(17) 



where p® is the asymptotic value of the background density. ip r denotes the real part of ^> g . Its imaginary part is v = v/c g is the 
soliton propagation speed measured in units of the sound speed c g for the GPE. It is given by 

c g = (Up° g /m)*. (18) 
In Eq. ( 17 1, Q < v < 1. Further, it can be shown that %jj r satisfies O 



(19) 



where = 1 — v 2 and w = (x — vt) /£, where the healing length £ is given by the well known expression £ = %/ J2mJJ p® — 
h/ \[2m c g . 



In analogy with the HGPE (see below Eq. ( 10 1), we define a dimensionless parameter A = c g /c , which is the sound velocity 



in the GPE measured in units of the zero-point velocity. This enables us to write £ in terms of A as 

£ = a/V2A. 



Hence w = \/2A(x — vt)/a — \/2Az. Using this in Eq. ( 19 1 and squaring both sides of this equation, we get 

dlp r 

dz 



(%) 2 = (l-^ 2 -z; 2 ) 2 



(20) 



(21) 
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FIG. 1: (color online) Condensate density profile p s (Eq. |12| ) of HGPE at half filling with V/t = 1/3, showing the behavior of the dark 
soliton as its propagation speed v — v/c s increases from to 1. Note that the soliton flattens out atti = 1. 



where z = Az. 



Now, from Eq. ( 17 1, the normalized condensate density p g = p g j p° for the GPE is 



Pg = i>f + y2 



(22) 



Differentiating Eq. d22l with respect to z, and using Eq. (21 1, it is easy to show that p g satisfies the following equation. 



( d jl? = {i-p g np 9 -v 2 ). 



(23) 
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The above equation describes the spatial evolution of the condensate density whose order parameter obeys the GPE, in a 
frame moving with velocity v. This equation is identical in form to Eq. ( 14 1 for the half-filled HGPE, and hence the condensate 
density p g for the GPE (which is equal to to the particle density) has the same functional form as the condensate density p s for 
the HGPE. We would like to emphasize that this mapping between the GPE and the HGPE is valid provided one considers only 
the condensate fraction of the HCB. Since the half-filling case for particle density in HGPE corresponds to quarter-filling for the 
condensate density, the corresponding GPE density for this mapping to be applicable is p a g = 1/4. 

The mapping between the GPE and the HGPE discussed above provides the relationship that the width of the soliton for the 
condensate density in the GPE and the HGPE are given, respectively, by T g = (2'yAy 1 and T s = (2-fA/Vl - 2A 2 ) -1 fl6l . 

We would like to note that the relationship between the GPE and the HGPE applies only to the condensate density profiles, 
as the solitons for the corresponding condensate wave functions in these two limiting cases will differ in their phase jumps. 



5. The BEC soliton and the magnetic soliton in the easy-plane ferromagnet 



An interesting consequence of the mapping between the GPE and the HGPE discussed in the previous section is the emergence 
of a relationship between two distinct physical systems, the GPE soliton in a BEC and a magnetic soliton that arises from the 
quantum easy-plane ferromagnet As we discuss below, this provides an interesting way to understand why there are in 
general two solitary waves in the bosonic density of the HCB system, and why in the special case of half-filling, the condensate 
part of the HCB has only one solitary wave and why it is dark. 

In magnetism, the study of strict solitons and solitary waves in various types of quasi-one dimensional classical spin systems 
has been a very active field ifTTl [T8l for over three decades. In the present work, the relationship between the Bose Hubbard 
model and a quantum ferromagnetic Heisenberg spin Hamiltonian ([3]) with an exchange anisotropy and a magnetic field will be 
exploited to relate the BEC soliton to the magnetic soliton. Note that the anisotropy is of the easy-plane type since g > 0. 

By writing the complex parameter r appearing in the spin coherent representation [15| as r = tan(6>/2) exp(i^), with 
< < ir and < <fi < 2tt, the coherent state average < S + > in the quantum spin- 1/2 system can be calculated in terms of 9 
and (f>. By identifying p s — \ < S + > | 2 we obtain the following relation between the condensate density p s and the classical 
spin variables: S x , S y ,S z , where 8 and <\> appear as the polar and azimuthal angles of a classical spin vector S. 



Ps = (S 2 X + Si)/ A = [(1/4) - 5 2 2 ]/4 (24) 

In addition, the particle density p, the chemical potential p and the particle-hole imbalance variable (1 — 2p°) in the BEC 
system described by HGPE are related to the spin variables of the magnetic system as follows. 



p = 1/2 -S z (25) 
fi = 2gp° = g(l - 25°) (26) 
g(l -2p°) = (g-p) = h z , (27) 

where h z is the magnetic field along the z-direction, as seen from the XXZ Hamiltonian ([3j. Of the above, the first equation 
is obtained by taking the expectation value of the corresponding hard-core boson operator to be spin coherent state average. The 
last two are obtained by using their asymptotic relationships in the expression for p (see below Eq. (j7]i). 

On taking spin-coherent state averages of the quantum spins on the lattice and going to the continuum, the evolution equation 
for the condensate order parameter also describes the evolution of classical spins, where the asymptotic density controls the 
external transverse field in the spin system. The solutions discussed here correspond to the half-filled case where this external 
field is tuned to zero. 



As a consequence of the relationship ( 25 i, a soliton solution of the total particle density p corresponds to a nonlinear excitation 



of the z-component of the spin. On the other hand, the soliton of the condensate density describes a nonlinear excitation of the 



in-plane spin, and relates to the square of the in-plane magnetization, due to the relationship (24i. 

As discussed in our earlier paper lfT2ll . the existence of the two solitons in HCB has its root in the fact that both particles 
and holes play equal roles in deciding the dynamics of the system. It is noteworthy that the XXZ Hamiltonian ([3) provides an 
alternative means to understand the existence of the two solitons. The spin Hamiltonian ([3]) contains terms with two competing 
effects, the easy plane anisotropy g that tends to align spins in the xy-plane and the the transverse field h z that favors a spin 
alignment along the z-axis. The ground state of the system consists of spins aligned on a cone that makes an angle 8q with the 
z-axis where 0q is determined by the background density of the atomic cloud given by po = 1/2 — S z = 1/2 — (1/2) cos 6q. 

Interestingly, in the half-filled case p° = 1 /2, the particle-hole imbalance Sq = (1 — 2p°) is zero, and h z vanishes. In this case, 
the ground state consists of spins in the xy-plane. Therefore, the solitary wave excitation resulting in a localized change in the 
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density makes the spins move out of the plane. In view of the fact that there is no clear preferred direction, disturbances above 
and below the easy plane are equally preferred. In the corresponding particle density in HCB, the dark and antidark solitons are 
therefore mirror images of each other, as seen from Eq. ( fTT) ). 



In contrast, as seen in Eq. (24 1, the condensate density of HCB is related to the square of the in plane magnetization, (which 
maps to the GPE density profile with p° = 1 /4) resulting in a single soliton both in the GPE and in the condensate fraction of 
the HCB. Such a soliton has to be a dark soliton, since the out of plane component of the spin decreases the total magnetization 
of the system. 

In summary, soliton characteristics in the BEC of bosons in weakly and strongly interacting regimes are expected to be quite 
distinct. In cold atom laboratories where the interaction between atoms can be varied, one expects significant changes in the 
dynamics as one tunes the scattering lengths of the bosonic atoms. Therefore, our result illustrating the similarities between the 
solitons in these two different regimes in the special half-filled case is important. Furthermore, by exploiting the relationship 
between the HCB and spin systems, the connection between the GPE and HGPE solitons paves the way for relating them to 
magnetic solitons in a ferromagnetic spin chain with easy-plane exchange anisotropy. Thus the HGPE soliton at half-filling 
and the above magnetic soliton can be regarded as two other incarnations of the GPE soliton. These studies open a new way 
to understand and interpret many important characteristics of solitons in these systems, including the fact that it provides an 
intuitive understanding of why the GPE soliton is dark. 

It is to be noted that it is the density in BEC that has a solitonic behavior. Thus by designing a quasi-one-dimensional optical 
lattice which incorporates the hard-core boson constraint of no double occupancy, and by loading them with bosonic atoms 
appropriately so as to simulate a half-filled lattice, it would be of interest to study soliton propagation of density in this strongly 
repulsive system, and investigate its relationship to the GPE dark soliton which is the hallmark of a weakly repulsive system. 
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